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1. Introduction 

Bose-Einstein condensate (BEC) with internal degrees of freedom has been much studied 
theoretically and experimentally. Bose-Einstein condensation of ultra-cold bosonic atoms has 
been realized under optical dipole traps. 1-3 This development stimulates a further research 
on physical properties of the condensate. Moreover, applications in atom optics 4 raise their 
potential, such as atom laser, atom interferometry and coherent atom transport. Internal 
degrees of freedom of atoms which are frozen under magnetic trap are liberated and play a 
significant role with wide possibilities in the condensate under optical traps. 

For BEC of dilute gases, the mean-field theory is valid. The time-evolution of the conden- 
sate is described by the Gross-Pitaevskii (GP) equation. It is known that the GP equation 
for single-component BEC is integrable in the one-dimensional (ID) uniform systems and 
has soliton solutions. In experiments, matter-wave dark and bright solitons are produced for 
single-component BECs. 5 ~ 8 

In general, we may suppose a system where each bosonic atom has the hyperfine spin 
of integer F. The condensate in the spin F state is described by the (2F + l)-component 
macroscopic wavefunction. Time-evolution of the multi-component condensate is derived from 
a generalized GP functional. The obtained multi-component GP equation should explain a 
variety of static and dynamical phenomena in the condensate. In particular, the existence of 
multi-component solitons is expected in the system. 

To analyze such nonlinear evolution equations exactly, we concentrate on the one- 
dimensional F = 1 spinor condensate. Recently, we have discovered a novel reduction to 
the integrable model. 9 ' 10 There, the atomic interaction is attractive and the spin-exchange 
interaction is ferromagnetic. In this paper, we consider the case of the repulsive atomic in- 
teraction and the anti-ferromagnetic spin-exchange interaction. We investigate in detail the 
one-soliton and two-soliton solutions and clarify their magnetic properties, which are absent 
in single-component solutions. We emphasize that the systematic study of multi-component 
soliton equations is rather new and has clarified interesting collision properties of solitons. 11,12 

The paper consists of the following. In § 2, we formulate the Gross-Pitaevskii equation 
for F = 1 spinor condensate, and present an integrable model with repulsive and anti- 
ferromagnetic interactions. The model is a reduction of a 2 x 2 matrix nonlinear Schrodinger 
equation. In § 3, solutions under constant boundary conditions are obtained by use of the 
inverse scattering method. Needless to say that the inverse scattering method for constant 
boundary conditions are much more involved than that for vanishing boundary conditions. 
The detail of the inverse scattering method is reported in a separated paper. 13 One-soliton 
states and two-soliton states are analyzed in detail respectively in § 4 and § 5. The last section 
is devoted to conclusion. The generic two-soliton solution is explicitly given in Appendix. 
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Bosonic atoms in the F = 1 hyperfine state are expressed as a three-component quantum 
field operator \l/ = (^i, ^q, ifr-i) whose equal-time commutation relations are 

[* a (x,t),ty(a/,t)] = 6 af ,6(x-a/), (1) 

for a,(3 = 1,0, —1. The interaction between atoms is supposed to be of short-range and have 
a form, 

V(xi - x 2 ) = S(xi - x 2 )(c + c 2 F 1 ■ F 2 ), (2) 

where Fi is the spin operator. Explicitly, the coupling constants are 

Co = 4:irh 2 (ao + 2a 2 )/3m, (3) 

c 2 = 4irh 2 (a 2 - a )/3m, (4) 

where a/ is the s-wave scattering length for the channel of total hyperfine spin / and m is 
the mass of the atom. / = 2 corresponds to parallel spin collision and / = to antiparallel 
spin collision. In this paper, we assume the system is quasi-one dimensional. The effective ID 
couplings are 

c = co/2a 2 j_, c 2 = c 2 /2a 2 j_, (5) 
where a± is the size of the transverse ground state. The second quantized Hamiltonian is 

As usual, the repeated indices mean the summation. Explicit form is the following, 

h 2 



(6) 



H = dx 



co + c 2 



d x ¥ a ■ d x V a + 
2m a 2 



+ *Li*Li*-i*-i + 2*J* (* f i*i + *ti*-i) 

+ (c - c 2 )^ t 1 ^ t _ 1 'i'_i^i + ^#5*5*o*o 



+ c 2 



(7) 
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In the mean-field theory of BEC, the dynamical variable $>(x,t) is the vacuum expectation 
value of the quantum field, 

*(x,t) = <*(M)> 

= (*i(M),$oOM),$-i(M)) T , (9) 

which is normalized to the total number of atoms Nt, 

dx&(x,t)-&(x,t) = N T . (10) 



The time-evolution of the spinor condensate wavefunction &(x,t) can be derived from the 
variational principle: 



3>i| 4 + |<I>_i| 4 



where the Gross-Pitaevskii energy functional is given by 
E GP = J 6x{^-d x ** a ■ d x $ a + 

+ 2|$ | 2 (|$i| 2 + |$-i| 2 )] + (co - c 2 )|$i| 2 |$-i| 2 

+ yl^ol 4 + P2(*i*Ii*o + $$ 2 $i$-i)}- (12) 

Substituting eq.(12) into eq.(ll), we get a set of nonlinear evolution equations for the spinor 
condensate wavef unctions: 

= - + (CO + C 2 )(|$l| 2 + |^o| 2 )^l 

+ (go - c 2 )|$_i| 2 $i + p^Ixfcjj, 

h 2 

ihdt^Q = - T^dl^o + (co + c 2 )(|$i| 2 + |$-i| 2 )$o 

+ co|$o| 2 $o + 2c 2 $o$i$_i, 
h 2 

\hd t <$>-i = - ^<9 2 $_i + (c + c 2 )(|$_i| 2 + |$ | 2 )$-i 

+ (c -c 2 )\<f> 1 \ 2 <f>- 1 +c 2 <f>* 1 <S>l (13) 

We call eqs.(13) multi-component Gross-Pitaevskii (GP) equation for F = 1 spinor Bose- 
Einstein condensate. 

Now, we consider a system with couplings cq = c 2 = c > 0. This corresponds to the situa- 
tion with repulsive inter-atomic interaction and anti-ferromagnetic spin-exchange interaction 
whose magnitudes are equal. In the dimensionless form: <& — ► (fa, \/20o, 0-i) T , where time 



and length are measured respectively in units of t = ha±/c and x = Hy a±/2mc, we rewrite 
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eqs.(13) as follows, 

idtfa = -dlfa + 2{\(j) l \ 2 + 2\U 2 )<Pi + tytrfl, 

id t <j>o = -& x <h + 2(|^i| 2 + \U 2 + l^i| 2 )0o + 2050i0_i, 

id t (f>-i = -d 2 x (t>-i + 2(|<£_i| 2 + 2|0 o | 2 )0-i + 2#>§. (14) 

These coupled equations are equivalent to a 2 x 2 matrix version of nonlinear Schrodinger 
equation (NLSE) with the self-defocusing nonlinearity: 

id t Q + 8 2 X Q - 2QQ^Q = O, 

with an identification, 

Q={ 4>1 <P °). 

Two remarks are in order. First, eq.(15) is a completely integrable system. Second, for unitary 
matrices U and V, if Q is a solution of eq.(15), then UQV is also a solution. 

3. Solitons of 2 X 2 Matrix NLSE with Nonvanishing Boundary Conditions 

In a separate paper, 13 we investigated the inverse scattering method (ISM) for the self- 
defocusing matrix NLSE with nonvanishing boundary conditions. It is shown that the standard 
form of the iV-soliton solution is expressed as 

N I n ie » y 

I + 2i(T^~I)S- 1 : . (17) 

\ U N e XN J 

Here / is the 2x2 unit matrix, and LTj {i = 1, • • • ,7V) are 2x2 Hermitian matrices, called 
the polarization matrices. S is a 2N x 2iV matrix defined by 



Q(x,t) =A e i ^*) 



e pi + e p 'j 

1 < i,j < N 



(18) 



where ifj G (0, 7r) and e Pj = sin ipj for j = 1, • • • , iV. The phase of the carrier wave is given 
by 

cj)(x, t) = kx- (k 2 + 2\l)t + S, 



and the coordinate function for the j-ih soliton is given by 

Xj(x, t) = — 2Aq sin <fj(x — 2(Aq cos (pj + k)t). 



(19) 



(20) 



Initial position of the j-th soliton is adjusted by scaling ILy. The analysis of the ISM shows 
that we should take detlLy = for all j. In that case, however, multi-soliton with the same 
spectral parameters turns to be one-soliton and is represented by a single polarization matrix 
with nonzero determinant (see Appendix). The expression (17) holds true even if this type of 
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degenerate multi-solitons are included. Therefore, we relax the condition so that detrL can 
have nonzero values. To repeat, Q(x,t) with arbitrary real detrL in general is the iV-soliton 
solution of eq.(15). 

The solution (17) has a form very close to the bright soliton solutions for the self- focusing 
matrix NLSE. However, there exists a significant difference: since the boundary conditions are 
nonvanishing and the nonlinearity is self-defocusing, the solitons are dark solitons, in general. 
The boundary condition in the limit x — > oo is fixed to be 

x — > oo. (21) 

The limit x — > — oo depends on whether detllj = or not, which will be discussed later for 
the one-soliton case and the two-soliton case. Other types of boundary conditions can be also 
realized from this standard solution by unitary transformations. 

For our setting of F = 1 spinor BEC, the explicit form of the iV-soliton solution is obtained 
through the identification (16) in eq.(17). Noting that Q in eq.(16) is a symmetric matrix, we 
naturally take LL to be real symmetric. An SU(2) transformation of a spinor wavefunction, 
<J>' = JJ<& with U = exp[i# ■ f], is equivalent to the following unitary transformation of Q: 

Q' = UQU T , U = exp[i0 • a/2], (22) 

where f = (P,f y ,F) are defined in eq.(8) and <r = (a x ,a y ,a z ) are the Pauli matrices. By 
this transformation, the boundary condition in the limit x — > oo can take any values except 
4>i = = 4>-i- For the exceptional case 4>i = = 4>-i, we may identify Q instead by 

«=L*!'V < 23 > 

V 9-1 4>o J 

Solitons in this case are just reduced to the one-component ones; (f>\ = = cf>-\ everywhere. 
That is, we realize one-component dark soliton solutions from the three-component solitons. 
Seen as a solution of the matrix NLSE, this corresponds to the degenerate one. For example, a 
twin-peak of the hole number density collapses to a single-peak. We shall discuss this situation 
in §4.2. 

The matrix NLSE has infinitely many conserved quantities due to its complete integrabil- 
ity. 11 First some conserved quantities are related to physical quantities as listed below: 

Total hole number: N T = J dx n(x,t), (24) 

n(x, t) = tv(QlQ ± ) - tr(Q ] Q). (25) 

Total spin: F T = J dx f(x,t), (26) 

f{x,t) =tr(Qt (X Q). (27) 

Total hole momentum: Pr = J dx p(x,t), (28) 
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p(x, t) = -ih[tv(QiQ ±jX ) - ti(Q*Q x )]. (29) 

Total hole energy: Et = j dx e(x,t), (30) 
e(x,t) = c[tr(Qi iX Q±, x + QiQ±QiQ±) 

-tr{QlQz + Q*QQ*Q)]- (31) 

Since there remain carrier waves at infinity, meaningful finite quantities are calculated by 
subtracting the background Q± = \im x ^± 00 Q(x). The bar denotes such a hole contribution 
calculated from the background. Physical quantities without the bar express particle prop- 
erties. Among local quantities, the spin density / = (f x ,f y ,f z ) is covariant under SU(2) 
transformation (22) and rotated as 

f' = Kf = ix{QW(rUQ). (32) 

If we choose U = exp[i7<7*/2] exp[i/3cr x /2] ex.p[iaa z /2] (a, (3, 7: Euler angles), the spin-rotation 
corresponds to 

U = K(a, (3, 7) = n z (-7)K x (-P)K z (-a), (33) 

where TZi(^) (i = x,y,z) is the matrix representing the rotation of an angle £ around the 
i-th axis. The other densities such as n(x,t), p(x,t) and e(x,t) are invariant under SU(2) 
transformation. The SU(2) symmetry of the system is attributed to the energy degeneracy 
for this spin rotation. The appearance of the spin density is the most remarkable feature for 
multi-component solitons. 

4. One-Soliton States 

In this section, we classify the one-soliton solution for our F = 1 spinor BEC. Setting 
TV = 1 in eq.(17), the standard form is expressed as 

Q(x,t) = A e i ^)S( X ;^n) 

= A i*(* )t ) (e 2 (*+f +i ^ det n + e^trll + 1)/ + 2ie* +i m 

e 2 (*+p) det n + e^+Ptrn + 1 ' 1 ' 

where we put e~ p = sin ip. Again, the phase of the carrier wave is given by 

4>{x, t) = kx - (k 2 + 2Aq)£ + 5, (35) 

and the coordinate of the envelop-soliton is given by 

x(x, t) = — 2Ao smip(x — 2(X cos (p + k)t). (36) 

In eq.(34), II is a real symmetric 2x2 matrix. In order to have a non-singular solution, it is 
necessary that det II > and trll > 0. 

There are 3 parameters and a matrix for a one-soliton: 

k : wave number of carrier wave, 
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Ao : amplitude of carrier wave, 

ip : the group velocity 2(Ao coscp + k) and 
the inverse of the size of soliton 2Ao sin 92, 

IT : contributions of the spin components. 

As is often the case for solitons, the amplitude determines the speed, and vice versa. The 
polarization matrix IT is a new ingredient for multi-component solitons and actually there is 
no counterpart in one-component solitons. The densities are explicitly calculated as follows, 

n(x, t) = 7 (e 2x+2p det II trll + 4e x+p det II + trll), (37) 

/(»,«) = 4> *y c w> «fc"'-i). <»> 

p(x, t) = h(k + Ao cos ip)n, (39) 

where det S = e 2 * +2 P det II + e^trl! + 1. 

We have two kinds of solutions corresponding to det II = and det II > 0, respectively. 
For detll = 0, the total spin is nonzero, while for det II > 0, the total spin is zero. Other 
than the total spin, the boundary conditions are also different for these cases. Therefore, we 
can distinguish one-solitons by checking the boundary values. We detail these two solutions 
in the next subsections. 
4-1 Ferromagnetic state 

Consider the case where det II = 0. We can write the symmetric matrix II as 

/ cos 2 § sin % cos f \ , N 

n = A 6 6 2 9' (40) 

\ sin § cos § sin z § / 

with a non-zero constant A. The boundary conditions of the standard form (34) are 

Qe"^ -► Ao/, x -> 00, (41) 
Qe~ i<p A UU T , (42) 

where 

U = e^/ 2 exp [-i6a y / 2] exp [iipa z /2] . (43) 

Comparing the left and right boundary values, we see that, not only the global phase but 
also the amplitude of each component are different. That is, the SU(2) rotated boundary 
conditions. It is remarkable that the profiles of the wavefunctions are in the shape of domain- 
walls. In contrast, the particle number density profiles have a dark soliton shape. See the 
upper panels of Fig.l. In particular, the spin is nonzero in total. Therefore, we call the soliton 
being in the ferromagnetic state. 
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Conserved quantities are integrated as follows. 

N T = 2\ sm(p, (44) 

- trllcr - rr 

F T = -N T —— = -N T (sm 6, 0, cos 9) T , (45) 
trli 

P T = N T h(k + X cos if) , (46) 
E T = f + V 

= N T c ((k + A cos ip) 2 + A 2 , - ^A 2 , sm 2 <p\ , (47) 
f = N T c (^(k + A cos iff -\ 2 Q + ^Xl sin 2 ^ , (48) 
V = N T c Uxl -^Xl sin 2 <pY (49) 



Changing the boundary conditions by global SU(2) transformation, we can vary the direction 
of the spin density as shown in eq.(32). 
4-2 Polar state 

Consider the case where detn > 0. In fact, the one-soliton with detll > is obtained as 
a multi-ferromagnetic-soliton solution where det 11^=0 and tpj are set the same value for all 
j. In spite of the multi-soliton content, it behaves as a one-soliton. The boundary conditions 
of the standard form (34) are 

Qe~^ A /, x -> 00, (50) 

Qe~^ A e 2iv /, x -00. (51) 

Unlike the ferromagnetic state, boundary values differ only by phase, i.e., U(l) transformation. 
One example of the density profiles is shown in lower panels of Fig.l. Observing the shape 
of the profiles, <j>±\ are dark solitons, i.e. hole-like envelopes, and fa is a bright soliton, i.e. a 
peak. The valleys of 4>±\ are of the same depth but at different positions and there is a peak 
of 4>o at the midpoint of them. The profile of the particle number density has a twin-valley. 
The profile of the spin density is dipole-like shape, i.e. an odd function of the coordinate, so 
that the total spin amounts to zero. For this reason, we call the soliton being in the polar 
state. By integration, conserved quantities are obtained as follows. 

N T = 4A sin <p, (52) 

F T = (0,0,0) T , (53) 

P T = N T h(k + A cos <p) , (54) 
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E T = f + V 

= N T c ((k + A cos ip) 2 + Aq - ^A 2 , sin 2 93^ , (55) 

f = N T c (^k + A cos ip) 2 - X 2 + ^X 2 sin 2 ^ , (56) 

V = N T c (2XI- I X 2 sin 2 ■ (57) 



The twin- valley shape merges when detll = (trIT) 2 /4, or equivalently II oc I. In this 
case, the diagonal elements of Q have the same values and the off-diagonal elements are zero 
everywhere. In other words, the soliton is reduced to the one-component one. In fact, the 
physical densities are all the same as those of the one-component NLSE soliton: 

n(x, t) = 2X 2 sinV sech 2 (( X + p + 0/2) , (58) 

f(x,t) = 0, (59) 

p(x, t) = h(k + Ao cos ip)n, (60) 

where we put e* = trII/2. Of course, the spin density is zero everywhere. 

One notices that if the left-half part of the polar soliton in the lower figures of Fig.l is 
pulled away to x = —00, then the remainder is the ferromagnetic soliton in the upper figures 
of Fig.l. This behavior is the same as that observed for the bright polar and ferromagnetic 
solitons in F = 1 spinor BEC with attractive and ferromagnetic interactions. 10 In fact, for 
fixed amplitude and group velocity, the total hole number, the total hole momentum and the 
total hole energy are just twice as those of the ferromagnetic soliton. It can be said that the 
polar state is a special case of two ferromagnetic states. 

5. Two-Soliton Collision 

In this section, we analyze the two-soliton solution in detail. The standard form of two- 
soliton solution is given by eq.(17) with N = 2. The explicit form is very complicated (see 
Appendix). 

The well-known and most remarkable property of solitons is the shape-preserving against 
mutual collisions. However, this is not always true for multi-component solitons in a naive 
sense. We look into the following three cases of two-soliton collisions: (i) polar-polar, (ii) 
ferro-polar, (iii) ferro-ferro. We concentrate on the case where soliton 1 moves rightward 
faster than soliton 2; soliton 1 will eventually pass over soliton 2. Each soliton moves with 
a constant speed except in the collision region. When another soliton is far apart, a soliton 
almost has the form of the one-soliton solution. To figure out the asymptotic form of soliton 
1 in t = ±00, fix the coordinate \i of the envelope of soliton 1 and then take the limit 
t — > ±00 44> e X2 — > or 00. Similarly, the asymptotic form of soliton 2 is obtained. 
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(a) (b) (c) 

Fig. 1. Snapshots of 1-soliton density profiles. The upper panels show plots for the ferromagnetic 
state with Ao = 1, k = 0, ip = 0.927, IT = ^ ^j. The lower panels show plots for the polar 

/ 8 3 88\ 

state with Ao = 1, k = 0, ip = 0.927, IT = I ' . (a) density of each component, \<j)i\ 2 (solid 

\3.88 2 / 

line), \(f>o\ 2 (chain line) and |</>-i| 2 (dotted line), (b) particle number density n, (c) spin densities, 
f x (solid line) and f z (dashed line). f y is zero everywhere for both cases. 



5.1 Polar state vs. Polar state 

Let us examine the collision of two solitons both in the polar state. We set detIL,- > 
for j = 1,2. Two solitons in each component collide each other in the same manner as the 
one-component two-soliton collision. We see that the only change due to the collision is a 
phase shift. Calculating the asymptotic forms from eq.(A-l), we have for t — > — oo, 

QT = A e^ +2i ^D (xi + s; pi, III) , (61) 
Q^ = A e i ^(x 2 ;P2,n 2 ), (62) 



and for t — > oo, 



Q~Qi n + Qt n , 

Qf n = A e i *D(xi;^i,n 1 ), (63) 
Qf = A e^ +2i ^!D ( X 2 + s; n 2 ) • (64) 
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The superscripts in and fin mean initial and final, respectively. Recall that D (%; (p, II) is 
defined in eq.(34). The position shift is determined by e s = S where 

„ _ sin 2 (5(^2-^1)) 



(65) 

sm 



The boundary conditions of the asymptotic forms are consistent for the two solitons. The 
results show the following: the soliton at the left hand is only changed by a phase shift and 
a position shift from the standard one-soliton solution in the presence of the soliton at the 
right hand which is in the standard form. The conserved quantities are calculated to be the 
sum of those of two asymptotic one-solitons. 
5.2 Ferromagnetic state vs. Polar state 

Let us see the collision between soliton 1 in the ferromagnetic state and soliton 2 in the 
polar state. We set det = and det IT 2 > 0. We write the polarization matrix II x for soliton 
1 as 

/ cos 2 % sin %■ cos % \ 
\ sm -y cos -j- sm -f J 
and accordingly we write (cf. eq.(43)) 

Wi = e ivi exp [-10^/2] exp [vpia z /2] . (67) 

Then, we calculate the asymptotic forms. The results are for t — ► —00, 

Q^QT + Qf, 

QT = Aoe^+^D (xi + s; ^ITJ , (68) 
Q^ = Aoe i *D( X2 ;^2,n 2 ), (69) 



and for t — > 00, 



Q^Q^ + Qf\ 



Qf n = A e i *D( X i;^i,n 1 ), (70) 
Qf = Aoe^iS) ( X2 ; n 2 ) U\. (71) 
Here, with H in eq.(65), we have introduced a deformed polarization matrix, 
~J Htrtnxna] sVatrK-ia^nxna] \ 

2 ^sV2 t r [n!(-i^) T n 2 ] tr [(-i^)n 1 (-iaf) T n 2 ] J ' 

In the presence of the soliton living at the right hand, the soliton living at the left hand has 
deformed boundary conditions. The deformed polarization matrix comes from the effects of 
the collision as well as the presence of another soliton at the right hand. The same is true for 
the ferro-ferro collision as will be seen later. 

Figure 2 shows the behavior of the ferro-polar collision. Notice that the time direction is 
downward. It is observed for soliton 2 (polar) that much of the hole amplitude is moved from 
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<fr±i to (j)Q due to the collision. For soliton 1 (ferro), no such a spin-mixing (or spin-transfer) 
is observed and its domain-wall shape is preserved. This is explained as follows. Soliton 1 has 
nonzero spin so it can rotate the local spin density of soliton 2; soliton 2 has zero spin in total 
so it cannot affect the local spin density of soliton 1. This kind of phenomenon was called a 
spin-switching in the attractive case. 9 ' 10 It is interesting to recognize that the spin density f x 
of soliton 2 shows up with a dipole-like shape after the collision (Fig.2.(d)). 



\4>i\\ 




(a) 



(b) 



(c) 



(d) 



Fig. 2. Density profiles of 2-soliton collision. Soliton 1 trailing right-down is ferromagnetic and soliton 
2 trailing left-down is polar. The parameters are fixed as Ao = 1, k — 0, ipi — ir — if2 = 1.12, IF = 

^ , n 2 = ^ . (a) density |</>i| 2 , (b) density |</>o| 2 , (c) particle number density n, (d) spin 

density f x . \(f>-i\ 2 behaves almost like |0i| 2 . f y and f z are zero everywhere. 



5.3 Ferromagnetic state vs. Ferromagnetic state 

Let us examine the collision of two solitons both in the ferromagnetic state. For j = 1, 2, 
set det IF, =0, and write 

IF, = exp [-Wj(T v /2] diag(Aj , 0) exp [-W j0 - y /2] , (73) 

and 

Uj = exp [-WjO-y /2] exp [i^-<7*/2] . (74) 
Then we have the asymptotic forms as follows. From eq.(A-l), we have for t —>■ — oo, 

Q^Q™ + Q™, 

QT = Aoe^aS) (xi! Vi&i) ^1, (75) 
Q^ n = Aoe i ^(x 2 ;¥'2,n 2 ), (76) 
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and for t — > oo, 

Qf n = A e i *S)(xi;yi,ni), (77) 

Q| n = Aoe^iS) ( X2 ; <P2, n 2 ) Wf, (78) 
where deformed polarization matrices are obtained as 

Uj = exp [-i6>V/2] diag(A, 0) exp [-10V72] , (79) 
for j = 1, 2. Here we have introduced new parameters, 

A = ~ cos 2 1(0 2 - 0x) + sin 2 1(0 2 - (80) 

tan ^=H-V2 tan 1(02-^). (81) 
We check the boundary conditions for the asymptotic forms. We have for t — ► — oo, 

Qfe -i * -»• A W 2 W 2 T , x ^ oo, (82) 

Qfe-W A WW T , x ^ -oo, (83) 

Q 2 n e"^ -»• A /, x ^ oo, (84) 

Q 2 n e~^ -»• A W 2 W 2 T , x -► -oo, (85) 



and for t — ► oo, 



where 



Q™e -i * -► A I, x ^ oo, (86) 

Qi n e-^ -► A WiWf, x -► -oo, (87) 

g 2 n e"^ -► A UiUf, x -> oo, (88) 

Q 2 n e"^ -► A WW T , x ^ -oo, (89) 

U = e i ^+v)/ 2 exp [-102*72] exp [w72] 

x exp [-i0V72] exp [i^ia z /2] (90) 
=e i(^ 1+ ^ 2 )/2 exp [_i0 1(7 2//2] exp [^1(^/2] 

x exp [-i0V72] exp [ic^72] . (91) 



An example of the profiles of two solitons before and after the collision is shown in Fig. 

3. The collision causes the spin precession. Let us see this phenomenon closely. Using eq.(32), 

we find that the spins of the two solitons have the following form, 

trTT™«- 

Fr = -\Fn—^r, |^ m |=2A sin^, (92) 
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Fig. 3. Snapshots before and after the collision of two solitons both in ferromagnetic state. The 
parameters are fixed as Ao = 1, k = 0, ipi = it — tp 2 = 1.12, IF = ^ j , II2 = ^ ^ . The 

collision takes place around t — 0. The upper (lower) panels depict the behaviors before (after) 
the collision at t = — 4 (t = 4). (a) density of each component, \(f>i\ 2 (solid line), |0 O | 2 (chain 
line) and |</>-i| 2 (dotted line), (b) particle number density n, (c) spin densities, f x (solid line), f y 
(chain line) and f z (dotted line). Soliton 1 and soliton 2 have spins F{ n = (—0.575, —1.19, 1.23), 
= (0,0,-1.8), Ff n = (-1.8,0,0), F 2 fin = (1.23,-1.19,-0.575). The rotation angle of spins 
around the total spin is w = —1.38. 



for i = 1, 2 and m = in, fin. We find that the polarization matrices are factorized as 



(93) 



where 



u m = e W2 exp [-i0 ja v/2] expliipja* /2] 



sm • 



for (i,j,m) = (1,2, in), (2,1, fin), and 



(94) 



u 



m _ / COS f 



sm y 



for (i,m) = (l,fin), (2, in). Explicitly, the spins are calculated as 



/ COS On 



'j cos ifj sin 9' + sin 9j cos 9' \ 



— sm ipj sm t) 
\ — sin 9j cos </?j sin 9' + cos 0j cos 9' y 



(95) 



(96) 
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for (i,j,m) = (1,2, in), (2,1, fin), and 

/ sin 6^ \ 

F l m = -\F™\ , (97) 
\ cos 9i J 

for (i,m) = (l,fin), (2, in). Also, for the total spin, we have 

\F T \ = 2A (sinVi + sin 2 y?2 + 2sin^i simp 2 cos6') 1/2 . (98) 

From eqs.(94) and (95), one can show the relation, 

< = -A (uf - (1 + e-'^E^ 2 )^ • <)« 3 j , (99) 

for (i, j, m, n) = (1, 2, fin, in), (2, 1, in, fin). Using the identity, 

tr[(ti£ • ui)« 2 • u\a] = - [Fx + F 2 + i(F 2 x F{)] , (100) 

for Fi = tr[«j • u\<t\ together with eqs.(92) and (93), we arrive at a formula for the rotation 
of spins, 

Ff J-{[(1 + A j7 ) + (1 - A j7 ) cos^]Fj n 
+ [(l + A J -,)-(l + A J -,)cosw]F, 



m 



F m x pin . 



-;- -2^m ' j ,2 * l }, (101) 



for (j, Z) = (1,2), (2,1) where 



\F T \ 



I p |2 _ I p |2 



The rotation angle uj around the total spin Ft is determined by 

sinw = — — ; — — — — - — - — , . (103) 

Recall that E is defined in eq.(65). 
6. Conclusion 

In this paper, we have studied soliton solutions of the integrable multi-component Gross- 
Pitaevskii equation for F = 1 spinor Bose-Einstein condensate (BEC) with repulsive and anti- 
ferromagnetic interactions. The solutions are directly obtained by identification with those of 
the self-defocusing matrix nonlinear Schrodinger equation (NLSE), which was solved by the 
inverse scattering method (ISM) in ref. 13. Since the interactions are repulsive, the solitons 
arise as dark solitons in general under nonvanishing boundary conditions. These dark solitons 
have similar properties compared to bright solitons in F = 1 spinor BEC with attractive and 
ferromagnetic interactions. 9,10 

We have shown that one-solitons are either in the ferromagnetic state or in the polar state. 
A soliton in the polar state is proved to be a special case of two solitons in the ferromagnetic 
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state moving with the same velocity. These two states, ferromagnetic or polar, are selected by 
choosing the boundary conditions. When one takes the left and right boundary amplitudes 
of wavefunctions differently, i.e., the SU(2) "rotated" boundary conditions, the one-soliton is 
in the ferromagnetic state. On the other hand, when the boundary conditions differ only by 
phase, we get the polar state. As is always the case for dark solitons, these changes in the left 
and right boundary conditions make solitons the topological objects (excitations), guarantee- 
ing the stabilities and the existence of finite spin amplitudes even in the anti-ferromagnetic 
interactions. Two-soliton collisions for every combination of the states have been examined. 
A soliton in the ferromagnetic state, which has nonzero total spin, causes spin rotation of 
another soliton by the collision. However, a soliton in the polar state do not rotate the spin of 
another soliton because it has zero spin in total. It can be said that, as "magnetic" carriers, 
solitons in the ferromagnetic state are operative while those in the polar state are passive. 

To detail the A-soliton collisions is straightforward. As is confirmed in the asymptotic 
forms of the two-soliton solution, when the other solitons are far apart, a soliton keeps the 
form of one-soliton with deformed boundary conditions and a deformed polarization matrix, 
affected by the solitons living at the right hand side of it. The explicit expression of the formula 
for multi-soliton collisions may be complicated and needs a further analysis. As investigated 
for the Manakov model in ref. 12, it is important to clarify the factorization property of the 
iV-soliton collision. 

Very recently, plane wave solutions and soliton solutions are obtained for general coupling 
constants cq and C2- 14 The stability of the solitons is numerically studied. 15 An interesting 
question is whether the multi-component Gross-Pitaevskii equation is integrable only at some 
special points as considered here. In particular, its generalization into the generic hyperfine 
spin F case is very interesting 16 and is left for a future problem. 



Appendix: Explicit Expression for Two-Soliton Solution 

We give the explicit form of 2-soliton solution of the matrix NLSE (15) as follows. 




(A-l) 



+ e pi+p2 trllitrll2 - sin 



-2 
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+ trIIie Xl+pl + detn 2 e 2(x2+p2) + trn 2 e X2+P2 + 1, 
and 23 = (QS^j) is a 2 x 2 matrix such that 
®n =H 2 e 2i(vi+V2) detni detn 2 e 2 ^ 1+pl ) +2 ^ 2+ ^) 



Full Paper 



+ S det ITi 
+ E det n 2 



e 2i ^ 1+ ^)(n 2 ) 11 + e 2i< ^ 1 (n 2 ) 

e 2i ^ 1+ ^)(ni) n +e 2i ^ 2 (ni) 



e 2(xi+Pi)+(X2+p2) _|_ e 2i<pi d e tnie 2( ^ Xl+ '' 1 ' ) 
P (xi+Pl)+2(X2+P2) 



+ {se a(,pi+,P2) (ni)ii(n 2 )ii + s(ni)22(n 2 ) 22 + e**' 1 (111)11(112)22 + e 2i ^(ni) 2 2(n 2 )ii 
- sin- 2 Q(^i + V2 )) e -(pi+p 2 )+K^i+^) [ ( ni)i 2 (n 2 ) 2 i + (ni) 2 i(n 2 )i 2 ] } e (»+p0+(x2+p 2 ) 

+ [e 2i ^(IIi)ii + (ni) 22 ] e x 1+ P! 

+ e 2 ^ 2 detn 2 e 2 ^ 2+ ^) + [e 2i ^ 2 (n 2 )ii + (n 2 ) 22 ] e* 2+ " 2 + 1, 



<Bi2 =2i{~e 2 ^ +2i ^ + ^ 2 detni(n 2 )i 2 e 2 » +X2 + ~e 2p2+iVl+2iV2 (U^u det U 2 e Xl+2x2 
+ [-S 1 / 2 e^ 2+i ^ +i ^ 2 (ni)i 2 trn 2 + S 1 / 2 e^ +i ^ + ^ 2 trni(n 2 )i 2 ] e Xl+X2 

+ e^(ni)i2e xi +e^ 2 (n 2 )i 2 e X2 }. 

We note that 5321 (resp. 2522) is given by replacing the indices of IL, with 1 <-> 2 in 25i2 (resp. 
®n). 

When we choose ipi = ip 2 = f for two ferromagnetic solitons where detlli = detll2 = 0, 
we have a one-soliton solution (34), 

Q(x,t)=A e i ^D( X ;^n), 

where Xi = X2 = X an d II = ITi + II2. This one-soliton is in the polar state if det II 7^ 0. 
Remark that det 111 = and det II2 = do not make det II = for II = IIi + LI2 in general. 
This kind of reduction gives multi-solitons in ferromagnetic and polar states. Then, we can 
consider (A-l) with det IT,- > and trllj > (j = 1,2) as the general form of 2-soliton 
solution. 
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